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' The heating of the solar corona is hkely to be due to reconnection of the highly com- 

\ plex magnetic field that threads throughout its volume. We have run a numerical 

I ' experiment of an elementary interaction between the magnetic field of two photo- 

\ spheric sources in an overlying field that represents a fundamental building block 

■ of the coronal heating process. The key to explaining where, how and how much 
\ energy is released during such an interaction is to calculate the resulting evolution 

T-H ■ of the magnetic skeleton. A skeleton is essentially the web of magnetic flux surfaces 

\ (called separatrix surfaces) that separate the coronal volume into topologically dis- 

■ tinct parts. For the first time the skeleton of the magnetic field in a 3D numerical 
\ MHD experiment is calculated and carefully analysed, as are the ways in which 

■ it bifurcates into different topologies. A change in topology normally changes the 
' number of magnetic reconnection sites. 

In our experiment, the magnetic field evolves through a total of six distinct 

\ topologies. Initially, no magnetic flux joins the two sources. Then a new type of 

■ bifurcation, called a global double-separator bifurcation, takes place: this bifurcation 
is likely to be one of the main ways in which new separators are created in the corona 

Q ■ (separators are field lines at which 3D reconnection takes place). This is the first of 

| H \ five bifurcations in which the skeleton becomes progressively more complex before 

^ ■ simplifying. Surprisingly, for such a simple initial state, at the peak of complexity 
there are five separators and eight flux domains present. 
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1. Introduction 

The solar corona, the upper atmosphere of the Sun, is known to be a low-density 
plasma with a temperature greater than 1 x 10® K, while the surface of the Sun - 
the photosphere - has a temperature of about 6000K. The solar atmosphere consists 
of active areas around sunspots, and regions of lower activity known as the quiet 
Sun. The quiet-Sun corona is, however, far from quiet and consists of magnetically 
open coronal holes, many types of coronal loop and dots of emission with a range of 
scales, such as microflares and X-ray bright points. Although the heating of X-ray 
bright points, and most probably microflares, has been shown to be produced by 
magnetic reconnection (Priest et al. 1994; Parnell at al. 1994a, b), the heating of 
coronal holes and coronal loops remains an unsolved mystery. 
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It has long been known that the heating mechanisms depend on the magnetic 
fields that are present over the surface and throughout the atmosphere of the Sun. 
For many years these mechanisms were generally divided into two categories: mag- 
netic waves and magnetic reconnection. Although low-frequency waves have now 
been observed in the corona (e.g., Nakariakov et al. 1999; De Moortel et al. 2002) 
and are important for coronal seismology (Roberts 1984; Dc Moortel 2005; Nakari- 
akov & Verwichte 2005), they do not carry enough energy to heat the corona and 
so reconnection is now regarded as the most likely mechanism. 

Various reconnection mechanisms have been proposed, including those that re- 
quire the presence of nulls (e.g., Lau & Finn 1990; Priest & Titov 1996; Longcope 
2001; Priest et al. 2003) and those that can occur in the absence of nulls (e.g., 
Parker 1972; Priest & Demoulin 1995; Birn et al. 1998). 

The quiet-Sun photosphere is covered by small magnetic fragments, which are 
driven across the surface by photospheric flows. New flux appears through a pro- 
cess called emergence which creates two opposite polarity fragments forming an 
ephemeral region (Harvey & Martin 1973). Flux is later removed through the pro- 
cess of cancellation when two opposite-polarity fragments collide and mutually lose 
flux (Martin et al. 1985; Harvey 1985). Cancellation and emergence lead to recon- 
nection that can power an X-ray bright point as explained by Priest et al. (1994), 
Parnell et al. (1994a, b) and von Rekowski et al. (2006a,b). However, the simple 
moving around of magnetic fragments in the photosphere also causes magnetic re- 
connection in the atmosphere (Longcope 1998; Galsgaard et al. 2000; Parnell & 
Galsgaard 2004; Galsgaard & Parnell 2005). The ceaseless movement of photo- 
spheric fragments in the quiet Sun drives a wealth of reconnection which can power 
X-ray bright points, microflares and nanoflares. Indeed, Close et al. (2004, 2005a) 
determined a conservative estimate that coronal magnetic flelds are recycled ten 
times faster (in just 1.4 hrs) than photospheric magnetic fields, implying that there 
is considerably more reconnection occurring throughout the solar atmosphere than 
is driven by emergence and cancellation alone. 

In this paper, we use the term fly-by to describe the passing of nearby fragments 
where flux is conserved. Using a 3D MHD code, Galsgaard et al. (2000), Parnell & 
Galsgaard (2004) and Galsgaard & Parnell (2005) calculated the reconnection rate 
and energy release of a simple set of fly- by models. They considered two discrete 
source fragments of opposite polarity, situated in an overlying fleld. The sources 
were driven anti-parallel to each other by advection flows in the photosphere at 
various angles to the overlying field. The fiux through the base was conserved 
throughout. The motions caused the flux lobe from each source to be pulled under 
the flux lobe from the other source. A twisted current sheet was formed between the 
two flux lobes as the sources became connected. Later, new current sheets started 
near the ends of the connected flux region and extended along it as the sources 
were disconnecting. It was suggested that the connection/linking of the sources is 
caused by separator reconnection (see also Longcope 2001; Priest & Titov 1996) 
in the twisted current sheet surrounding this separator, and that the sources later 
disconnected through a process that they called separatrix-surface reconnection. 
Furthermore, the rates of reconnection were found to be fast: up to 58% of the 
instantaneous, perfect reconnection rate of the potential evolution. Although, the 
rates of reconnection were found to be almost a factor of two slower in the discon- 
necting phase than in the connecting phase. 
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Our primary aim is to understand much better where and how reconnection is 
occurring in the solar corona by first understanding in detail just how it is occurring 
in the fly-by numerical experiments of Galsgaard and ParncU. To shed light on this, 
we must first answer the following two fundamental questions in solar MHD: 

1. What is the topological structure of the magnetic field (i.e. where are the 
possible reconnection sites) as separate 3D magnetic flux systems interact in the 
solar atmosphere, and 

2. What kind of bifurcation leads to a change of topology (i.e. how and why do 
the number of reconnection sites change)? 

As a background, we need flrst to say a little about magnetic topologies and about 
3D magnetic reconnection. 

First of all, the topological structure of a complex magnetic field is best described 
in terms of its magnetic skeleton defined by e.g. Bungey et al. (1996) and Priest et 
al. (1997) on the basis of earlier work and definitions by e.g. Greene (1988), Lau & 
Finn (1990) and Priest & Titov (1996). Such a skeleton includes: 

1. sources, where the magnetic field enters or leaves the region being studied, 

2. null points where the magnetic field strength vanishes, 

3. separatrix surfaces, which arc made up of field lines that extend to or from a 
null and form the boundary between topologically distinct domains (also known as 
E or fan surfaces), 

4. spines, isolated field lines extending from a null back to a sink, if the separatrix 
surface from the null goes to a source, or vice versa (also known as 7 lines), 

5. separators, which are the intersections of two separatrix surfaces. A separator 
connects two null points and represents the dividing line between four topologically 
distinct domains (also known as A-B lines). 

Next we must define what a fiux domain is. A flux domain is a simply-connected 
volume of magnetic field lines that share the same start and end sources or source 
pair and is bounded by a separatrix surface (Longcope & Klapper 2002; Beveridge 
& Longcope 2005 and Parnell 2007). Thus, in such a domain it is possible to con- 
tinuously deform the field lines from one to another. This definition is related to 
the key feature that is generally used to find most flux domains, namely the con- 
nectivity of field lines, but it goes much further. If there were two separate paths 
linking the same source pair, each surrounded by their own separatrix surfaces, and 
if each of these volumes were simply connected then they would be classed as flux 
domains of a multiply- connected source pair. The main point here is that knowing 
the number of connected source pairs is not the same as knowing the topological 
structure of the field. 

In general, it is assumed that most source pairs have multiplicity one. Although 
as soon as the magnetic field becomes complex, as it is on the Sun, this is not 
necessarily the case. Source pairs with multiplicity two were found by Close et 
al. (2005b) and Barnes et al. (2005) in potential magnetic field extrapolations of 
two different magnetograms. In the present experiment, we find source pairs with 
multiplicity up to three. 

The simple potential topologies associated with three or four flux sources and the 
bifurcations between the topologies have been studied by Brown & Priest (1999), 
Beveridge et al. (2003) and Maclean et al. (2006). In particular, they discovered local 
bifurcations, in which the topology changes when pairs of null points are created 
or destroyed, and global bifurcations, in which the number of null points is held 
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fixed during a change of topology. A comprehensive review of such topologies has 
recently been given by Longcope (2005). Of particular interest to our numerical 
experiment is the Beveridge-Longcope equation (Beveridge & Longcope 2005 and 
Parnell 2007), namely, 

V = ^nDn = X + S-N^-l, (1.1) 

n 

which relates the number of separators (X), sources {S) and coronal nulls (A/^) to 
the number of domains {V = X]„»^-D„), where £>„ denotes the number of source 
pairs with multiplicity n. 

In this paper, we compare om- numerical MHD results with those for potential 
fields and discover that MHD evolution leads to a much richer and more complex 
set of topologies than evolution through a scries of equi-potential states. 

The paper is divided into the following parts. The next section, §2, is a descrip- 
tion of the potential and MHD models, while §3 describes the evolution of their 
skeletons. §4 and §5 highlight the nature of the bifurcations and the development of 
multiply-connected domains. Finally, in the last section, we discuss the implications 
of the results. 

2. Numerical Model 

The numerical model used here has the same set-up as in Galsgaard et al. (2000), 
Parnell & Galsgaard (2004) and Galsgaard & Parnell (2005). We give a summary 
here for completeness. 

The setup comprises a Cartesian box of 128 x 128 x 65 grid points with a scale 
of 1 X 1 X 1/4. The normal components of the initial magnetic field are imposed 
on the base [z = 0). They consist of a positive source, of radius Tq = 0.065 and 
maximum magnetic field strength Bq = 0.85, placed at (|,|,0) and a negative 
source with the same radius and maximum magnetic field strength, but opposite 
polarity, placed at (|, |,0). These sources have a cosine profile of the form = 
Bq (1 -|- cos(7rr/ro)) /2, where r < Tq is the distance from the centre of the source. 
Elsewhere on the base, the normal component of magnetic field is taken to be zero. 
Otherwise, the box is closed at the top and bottom and periodic in x and y. An 
overlying field, Bover = 0.12Boy is added to disconnect the two sources in the 
initial state. For comparison we consider a potential model and a dynamical MHD 
model which are described below. 

(a) Potential Model 

The potential magnetic field (Bp), by definition, satisfies V x Bp = and 
V.Bp = 0, so that it may be written as Bp = V^, where (j) satisfies 

V^4) = {). (2.1) 

Using a multigrid method with a Gauss-Seidel smoother and the above boundary 
conditions imposed on a unit box, we can solve equation 2.1 for (f) at each time 
frame, from which we evaluate Bp. The taller numerical box of height 1 for the 
potential field in comparison to the 1/4 for the MHD model makes essentially no 
diflEerence to the results as all the magnetic field at height 1/4 and above is simply 
the horizontal overlying field. 
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(b) Dynamic MHD Model 

The dynamic MHD approach uses the non-dimensionahsed, non-ideal MHD 
equations in the form: 



dt 



-V. (puu + f ) - Vp + J X B, 



— = - V.(eu) - pV.U — + Qvisc + Qjoule, 

C't ^cool 

_ ^ dp ^ , , 

E = - (u X B) + TyJ, J = V X B, 

P = (7-l)e, T=P, 

P 

where B, u, e, J, E, p, p, T, Tq, t, icooi, Q joule? ^visc cind T arc the magnetic 
field, velocity, internal energy, current, electric field, pressure, density, temperature, 
equilibrium temperature, time, exponential cooling time, magnetic resistivity, joule 
dissipation, viscous dissipation and viscous stress tensor, respectively. An ideal gas 
with the ratio of specific heats 7 = 5/3 is assumed. The equations have been non- 
dimensionalised such that the permeability of free space (/Uq) and the ratio of the gas 
constant over the mean molecular weight {R/fi) both disappear from the equations. 

Initially, we assume an equilibrium with a potential magnetic field in an atmo- 
sphere with uniform density and pressure of 1/2 and 1/6, respectively. The magnetic 
field is calculated using the method described in §2a using the two sources defined 
above in their initial state. 

These sources are advected rigidly by fiows in strips that are just wider than 
their corresponding sources, with the velocities ramped down along the edges of 
these strips to zero. The flows are accelerated from rest up to a constant speed of 
0.02 of the peak Alfven speed in approximately 0.6 of an Alfven time. The positive 
source is driven along the base in the x direction, and the negative source driven 
in the — x direction. 

The equations are solved using a third-order predictor-corrector method with 
a sixth-order method used to calculate the spatial derivatives. As staggered grids 
are applied, a fifth-order method is used for spatial interpolation. Viscosity and 
magnetic resistivity are handled using a fourth-order method combined with a dis- 
continuous capture mechanism to provide the highest possible spatial resolution 
for the given numerical resolution. Diffusion, and therefore magnetic reconnection, 
is anticipated to take place in current sheets that have a very narrow, fine-scale 
structure. In codes using uniform or anomalous current dependent resistivity (i.e. 
rj oc |j|) the resistivity (rj) has a fc^ dependence, where k is the wave number. 
This means that significant diffusion can occur even in well resolved current struc- 
tures. The advantage of the hyper-resistivity used here is that a higher k depen- 
dence (fc^) is achieved implying that the dissipation increases very rapidly as the 
wavelength decreases localising the diffusion more efficiently to short length-scales. 
This is equivalent to having finer-scale current sheets, but with the same numer- 
ical grid resolution. To achieve this the resistive and viscous dissipation are de- 
pendent on three factors, namely, (i) compression regions (for the resistivity this 
compression is only relevant when it acts perpendicular to the magnetic field vec- 
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tor), (ii) advection speed (plasma velocity) and (iii) the typical propagation speed 
of information (fast-mode/sound speed). More details of how these factors influ- 
ence the viscosity and resistivity and a description of the code can be found at 
http://www.astro.ku.dk/~kg/. 

Further particulars of the experiment can be found in Parnell & Galsgaard 
(2004). All the times in this paper are given in terms of the initial Alfven time of 
the experiment. 

3. Evolution of Magnetic Skeleton 

The corona has often been assumed to be near potential and many theoretical treat- 
ments have, for simplicity, assumed potential fields, the minimum-energy magnetic 
field for any normal field component prescribed. Potential fields may evolve through 
a series of equi-potential states by using perfect and instantaneous reconnection. It 
is useful, therefore, to compare the differences between the skeleton of our dynamic 
MHD magnetic fields with those that would arise from a potential field. 

For this, we classify the four difii'erent types of flux domain that are found in our 
experiment by way of their flux connectivity (or source pairs) . Flux in the overlying 
domain is not connected to either source (i.e. it connects one side boundary to 
another side boundary). The positive open{negative open) domain contains flux 
that connects the positive (negative) source to a side boimdary. Flux in the closed 
domain connects the positive source to the negative source. The boundaries of 
these flux domains deflne the skeleton. We discuss the evolution of skeletons of the 
potential and dynamic MHD models in turn. 

(a) Potential Skeleton 

The potential field is calculated using the method described in §2a for each time 
step, from which we determine the magnetic skeleton. We visualise the skeleton in 
3D (figure 1, col 1). The footprint of the skeleton represents the intersection of 
the 3D skeleton with the photosphere (figure 1, col 2). All null points found are 
photosplieric and the footprint is calculated by tracing the separatrix field lines 
that lie in the photospheric source plane from these null points. The spines are not 
shown. A cross-section of the skeleton at y = 0.5, which lies midway between both 
of the sources and the two null points, is also determined (figure 1, col 2). Both 
separatrix surfaces and all of the separators intersect this plane. 

By considering the evolution of the field, we group contiguous time frames into 
phases with the same configuration of flux domains. Four phases are found (see 
table 1): PI - initial open, P2 - closing, P3 - reopening and P4 - final open. We now 
describe these four phases in turn. 

Phase PI (0.00 < t < 1A7): Initial Open 

The first phase has three flux domains: overlying, positive open and negative 
open, where the positive source lies to the right of the negative source (flgure 1, 
row 1). There are no separators and no reconnection occurs. Thus all the flux is 
open, as seen in flgure 2. This phase lasts while the sources and their flux domains 
are pulled towards each other, but ends as soon as the domains touch. 
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Phose PI 




'Transition P'2/P3: 




Figure 1. Snapshots of the skeletons for the potential model at t = 0.0, 4.1, 8.14, 10.5, 19.9 
Alfven times. Column 1: 3D view, with positive separatrix surface (grey lines), nega- 
tive separatrix surface (black lines) and separators (white lines). Column 2: Footprint 
(0 < a; < 1, < J/ < 1) at z = 0, showing sources (hollow black circles), null points (posi- 
tive - V s-nd negative - A) and intersections of separatrix surfaces with the photosphere 
(solid lines). The dashed lines show the separators above the source plane for comparison. 
Column 3: Cross-section (0<a;<l,0<2< 1/5) at y = 0.5 showing the plane's intersec- 
tion with the separatrix surfaces from the positive null (thin line) and the negative null 
(thick line) and with the separator (o). 
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Table 1. Characteristics of the potential phases 

(The components of each phase for the potential model, namely, the number of separators 
(X), the multiplicity of each source pair and the total number of flux domains (©).) 
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Figure 2. Percentage of open (dashed), closed (solid) and reopened (dotted) flux from the 
negative source versus the time if) measured in Alfven times for the potential model (from 
Parnell & Galsgaard 2004). Vertical lines represent the transition times between phases. 



Phase P2 (1A7 < t < S.IA): Closing 

A global separator bifurcation occurs in which a separator rises from the source 

plane. This separator runs the ridge of a new domain that is made up of new recon- 
nected closed flux, (figure 1, row 2). Reconnection creates closed flux and destroys 
open flux throughout this phase (figure 2). The sources continue to move towards 
one another during this stage, until they reach their point of closest approach and 
are aligned with the overlying field. At this point the two sources are completely 
connected, the closed flux reaches its maximum (figure 2) and a new phase must 
start. 

Phase PS (8.14 < t < 15.49;.- Reopening 

As the separatrix surfaces of both nulls fully coincide, a global separatrix bifur- 
cation occurs in which an infinite number of separators completely encompasses the 
closed flux domain isolating it. This is a new type of bifurcation (figure 1, row 3). As 
soon as the sources move out of alignment with the overlying field the flux reopens 
through reconnection and a new single separator state is formed. The continued 
advection of the sources creates a positive open domain on the left and a negative 
open domain on the right with a closed domain in between (flgure 1, row 4). These 
domains are all simply connected and surrounded by overlying fleld. A separator 
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runs the ridge of the closed domain, as in phase P2. It descends towards the base as 
closed flux (with overlying flux) is converted back into positive and negative open 
flux (sec figure 2). This is the reverse of phase P2 when the closed flux was formed. 
This phase lasts until the separator has reached the base 

Phase P4 ("15.49 < t < 26.54): Final Open 

In this phase (figure 1, row 5) a global separator bifurcation destroys the sep- 
arator when it reaches the source plane, leaving just three simply connected flux 
domains: overlying, positive open and negative open. These flux domains are equiv- 
alent to those in phase PI, except that, as the sources have now passed each other, 
they are now moving apart instead of together and the negative source now lies to 
the right of the positive one. There is no reconnection in this phase (see figure 2). 

(6) Evolution of MHD skeleton 

The dynamic MHD experiment is executed using the method described in §2b. 
From this, the skeleton for each time frame is deduced. As in the potential model, 
we visualise the skeletons in 3D and also in 2D using both the footprint at 2; = and 
a cross-section at y = 0.5 (see figure 3). A filled contour plot of current intensity is 
added to the cross-section. 

As before, we group contiguous time frames into phases of the same topology. 
These phases are: Dl - initial open, D2 - double-separator hybrid, D3 - single- 
separator closing, D4 - quintuple-separator hybrid, D5 - triple-separator hybrid and 
D6 - single-separator reopening (table 2). We use a naming convention based on the 
number of separators and the reconnection process occurring, e.g. closing, hybrid 
(closing and reopening), reopening (see Parnell et al. 2007 for physical explanation). 
The phases are described in turn. 

Phase Dl (0.0 < t < 5.49): Initial Open 

Since the magnetic field is assumed to be potential at t = 0, phase Dl (figure 3, 
row 1) is topologically equivalent to phase PI. Hence, there are no separators and 
three simply connected flux domains, namely: overlying, positive open and negative 
open. The positive source (and its corresponding flux domain) lies to the right of 
the negative source (and its flux domain). 

As the positive open and negative open flux domains arc advectcd towards one 
another the two flux domains push up against each other, save for a small gap below 
(figure 3, row 1, col 2). At the time this phase ends, about half of the closed flux 
in the potential experiment has been reconnected (figure 2). Thus, since, up to this 
point, there has been no reconnection in the dynamic experiment, it is natural that 
a thin narrow current sheet has developed by the end of this phase. 

Phase D2 ("5.49 < t < 9.68J; Double- Separator Hybrid 

The two open flux domains break through each other marking the start of this 
new phase (flgure 3, row 2). Two separators and a closed flux domain are created by 

way of a global double-separator bifurcation (sec §4 for further details) . The closed 
flux domain is purely coronal and thus cannot be seen in the footprint, though it is 
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Table 2. Characteristics of the dynamical MHD phases 

(The components of each phase for the dynamic MHD experiment based on the number of 
separators [X) and the multiplicity of each source pair. The number in brackets denotes 
the numbers of domains that are coronal for each source pair. The total number of flux 
domains (2?) is also calculated.) 

Phase Seps. Source Pairs - Flux Domains Total 
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Figure 4. Percentage of open (dashed), closed (solid) and reopened (dotted) flux from 
the negative source versus the time {t) given in Alfven time units for the dynamic MHD 
experiment (from Parnell & Galsgaard 2004). Vertical lines denote the transition times 
between phases. 



clearly visible in the cross-section. This coronal domain pierces the overlying region 
creating a flxix region of multiplicity two. The field lines in the overlying flux region 
are no longer simply connected, but form two separate flux domains (table 2). 
The spines marked in the footprint denote the lines along which the positive and 
negative open flux domains fold back and touch the base, thus forming a separatrix 
surface boundary between the two overlying flux domains (see §5 for a detailed 
explanation) . 

A current sheet extends along the upper separator. This is visible as a tall thin 
current sheet in the cross-section (figure 3, row 2, col 3). Here, it is found that 
separator reconnection is creating closed magnetic field. However, around the lower 
separator reconnection is reopening the field and thus decreasing the flux in the 
trapped overlying flux domain. The shrinking of the small trapped overlying domain 
brings the lower separator to the photosphere signalling a change of topology. 

By the end of this phase 15% of the flux in the dynamic experiment is closed, 
whereas in the potential model all the flux has been closed and some 25% has 
reopened again. 
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Phase D3 <t< 16.96/- Single- Separator Closing 

A global separator bifurcation at the base destroys the lower separator and 
reduces the overlying field to a simply connected domain once again. Thus, phase 
D3 has just one separator and four flux domains, one of each connectivity (figure 3, 
row 3). From figure 4 it can be seen that the reduction from two to one separator 
has no effect upon the rate of closing. This is because the reconnection at the lower 
separator was weak. By the end of this phase, just 68% of the fiux is closed in this 
experiment, whereas all the flux in the potential model has both closed and then 
reopened, and the final phase P4 has begun. 

The separator in this experiment continues to create closed flux, which enlarges 
the closed flux domain until its encompassing separatrix surfaces expand out to 
touch the surfaces separating the open and overlying flux domains in two positions 
above the source plane (figure 3, row 3, col 3). 

Phase D4 (^16. 96 <t< 21.14): Quintuple- Separator Hybrid 

Here, two new separator pairs and four new coronal domains are created (fig- 
ure 3, row 4). The separator pairs are created through two global double-separator 
bifurcations. The bifurcations allow the positive (negative) open source pair to split 
into two through the creation of a new negative (positive) reopened flux domain 
generating the four new coronal domains. These six open flux domains can be seen 
in turn in the cross-section y = 0.5 forming an arc around the closed domain (fig- 
ure 3, row 4, col 3). From left to right these domains are: negative open (original), 
positive open (new re-opened), negative open (original), positive open (original), 
negative open (new reopened) and positive open (original). Hence, there are two 
source pairs, each of multiplicity three. Separating these flux domains are five sep- 
arators which all lie on the same pair of separatrix surfaces, and hence all connect 
the same two null points. 

The five separators represent five plausible sites for reconnection, but only three 
lie in high current regions (figure 3 row 4, col 3); thus significant reconnection is 
expected only at these sites. We find that reconnection at the central separator 
is closing the field whilst at the two separators on either side it is reopening the 
field. This behaviour is very different from the potential model where closing and 
reopening never happen concurrently. 

Reconnection at the lower separators converts the flux in the two small trian- 
gular open domains into overlying and closed flux, causing these separators to fall 
to the photosphere and the topology to change. 

Phase D5 (21.14 < t < 40.6): Triple- Separator Hybrid 

A pair of global separator bifurcations destroy the two lower separators, leaving 
three separators, and reduces the multiplicity of the positive and negative open 
source pairs to two. Thus, here, there are two source pairs of multiplicity one, 
namely: overlying and closed, and two of multiplicity two, namely: positive open 
and negative open. From figure 4 it is clear that there is still reconnection occurring 
at all the separators; the top separator is still closing the field, while the other two 
separators are reopening the field. 
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The driver is switched off at f = 25. Shortly after this we see that the symmetry 
in the system starts to be lost with the original negative fliix being used up much 
quicker than the original positive flux. Once the original negative domain is emptied 
of flux this phase ends and a new one beings. 

Phase D6 ^40. 6 <t < 56.4^; Single- Separator Reopening 

A global double-separator bifurcation again heralds the start of this phase. How- 
ever, here, wc do not sec the creation of a separator pair, but rather the loss of the 
left-hand and middle separators. This then leaves four flux domains, one of each 
connectivity (figure 3, row 6). This phase is topologically equivalent to phase P3 
of the potential experiment, although here the two scparatrix surfaces are rather 
distorted. Gradually symmetry is restored to the system and what was the right- 
hand separator now becomes a new central separator. The closed flux continues to 
be reopened along this separator until almost all the closed flux is used up. 

Just before the positive and negative open flux domains completely separate 
from one another the experiment ends. However, it is clear that if the experiment 
had been run for just a short time longer then a final phase D7 would have been 
entered, in which there would be three simply connected domains and no separators. 
This stage would have been topologically equivalent to phase Dl, but with the 
positive source and flux domain on the left rather than the right. It would be 
equivalent to phase P4 from the potential experiment. 

4. Global Bifurcations 

In an individual topological phase the number of reconnection sites remains con- 
stant, although their locations may change. A change in topology normally leads 
to an increase or decrease in the number of reconnection sites and occurs by way 
of a local or global bifurcation. During our experiment three global bifurcations 
arc found, two of which arc now. The global separator bifurcation, in which a single 
separator is formed above the source plane and a mirror separator is formed below, 
is well known (e.g.. Brown & Priest 1999; Beveridge et al. 2003 and Maclean et al. 
2006). During the potential evolution this bifurcation occurs between phases PI and 
P2 where it creates a new separator and it also occurs between phases P3 and P4 
where is destroys a separator. The bifurcation between phases P2 and P3, however, 
is new. It is a global bifurcation since the magnetic field in one phase cannot be 
deformed into another and it does not change the number of nulls or separators. We 
called it a global separatrix bifurcation. It changes the handedness or nature of the 
separatrix surfaces: before this bifurcation the positive fiux domain is to the right 
of the negative one, but afterward the order of the flux domains is swapped. Note, 
that this bifurcation is non-generic since it exists only when the two sources have 
exactly the same magnitude of flux: for example when the magnitudes arc slightly 
different it will be replaced by two global separator bifurcations, one to an enclosed 
topology and one away from it. 

In the dynamical MHD experiment, wc sec another new type of global bifur- 
cation in which a pair of separators that lie above the source plane is created or 
destroyed. We call it a global double- separator bifurcation. The separators are cre- 
ated by the bulging of the two separatrix surfaces through each other. This naturally 
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creates two curves of intersection between the two surfaces after they have touched 
and, hence, a pair of separators Unking the same pair of nulls. The reverse leads to 
the destruction of a separator pair. Wc suggest that this is one of the most natu- 
ral ways in which separators can be formed or destroyed in the solar atmosphere. 
Other possible ways are by a global separator bifurcation, a global spine-fan bi- 
furcation (Maclean ct al. 2005) and less commonly by local separator and local 
double-separator bifurcations. 

It is interesting to note that the reconnection that occurs at these pairs of 
separators each creates and destroys different types of flux. For instance, in phase 
D2, the upper separator reconnects positive and negative open flux to form closed 
and overlying flux. However, the lower separator reconnects closed flux with the 
overlying flux trapped below this domain to produce more positive and negative 
open flux. And similar opposite reconnection processes occur at the two pairs of new 
separators in phase D4. This means that for the new coronal domain bounded by 
the pair of separators to grow, the rates of reconnection along the separators must 
be different with the one filling the coronal domain dominant. It is natural that 
for a global double-separator bifurcation to destroy two separators the dominant 
reconnection process must be that which destroys the flux in the coronal domain 
bounded by the separator pair. 

The reverse reconnection along the new pair of separators is in contrast to 
the behaviour of the separator and mirror separator formed in a global separator 
bifurcation. Here, as one would expect from the name, the reconnection along the 
mirror separator, below the source plane, mirrors the behaviour along the separator 
above. That is their reconnection leads to the filling the domain they bound which 
straddles the source plane. Hence, global separator and global double-separator 
bifurcations show quite distinct behaviours. 

5. Multiply-Connected Domains and Magnetic Topology 

Formulae 

One of the most striking differences between the potential and the dynamical MHD 
evolution of the magnetic topology is the appearance of multiply-connected source 
pairs. These arc a natural consequence of the global double-separator bifurcation 
which pierces a hole in a domain, allowing the insertion of a newly created domain, 
and leading to distinct paths connecting the same pair of sources. Each flux domain 
of a multiply-connected source pair contains field lines that connect the same pair 
of sources and forms a simply-connected volume that is bounded by a separatrix 
surface. 



The global double-separator bifurcation occurs first during the dynamic MHD evo- 
lution at the start of phase D2 where it traps part of the overlying field whilst 
creating a new closed domain. Thus, the overlying field region is no longer a simply 
connected volume, but has two distinct paths, as discussed earlier. Prom the cross- 
section in the y ~ 0.5 plane (figure 3, row 2, col 3) it is clear that in this plane the 
entrapped portion of overlying field is bounded on the left by part of the separatrix 
surface from the negative null and on the right by part of the surface from the 
positive null; above, it is edged by a separator. However, what about in the y = 0.1 
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Figure 5. Cross-section at j/ = 0.1 of the magnetic field from the dynamical experiment at 
t = 7.6 Alfven times, i.e. during phase D2, showing the positive separatrix surface. Note, 
a;2:-plane only extends from (0.4,0.0) at the lower left to (0.55,0.005) at the upper right of 
the region. 

plane? Here, clearly, the separatrix surface from the negative null does not cross this 
plane, but its spine, which cannot divide topologically distinct regions as it is only 
a line, does. Focusing in on this region (figure 5) reveals that the positive separatrix 
surface is deformed and folds over such that it kisses the base at the spine from the 
negative null. Thus, the region of trapped overlying field is indeed bounded along 
its entire length by separatrix surfaces: between y — and the negative null it is 
enclosed by the positive null's separatrix surface alone; between the two nulls it is 
bounded by both the positive and negative nulls' separatrix surfaces; and between 
the positive null and y = 1 it is bounded by the negative null's separatrix surface 
alone. 

(a) Beveridge-Longcope Equation 

Since the global double-separator bifurcation does not occur during the po- 
tential evolution of the magnetic field, all the potential source pairs have multi- 
plicity one. Thus, for all the potential phases, the sum over all the source pairs 
V = J2n=i iT-Dn = Di. From table 1 we see that for every generic phase the differ- 
ence between T> and the number of separators (X) is 3. This is a consequence of 
the Beveridge-Longcope equation (equation 1.1), with four sources, including two 
for the overlying field (5 = 4), and no coronal nulls (iVg = 0). 

In the first phase of the dynamic experiment, Dl, has the same topology as phase 
PI, hence it too must satisfy the Beveridge-Longcope equation. From table 2, phase 
D2 has 3 source pairs with multiplicity one and one source pair with multiplicity 
two, so here V = ^n=i'"'^n = 1x3-1-2x1 = 5. Thus, from equation 1.1, 
T) ~ X — S — Nc — 1 = 3, since there are two separators, X = 2. As in the potential 
case, of course, 8 = 4 and A^c = 0, since the number of sources has not changed 
and again no coronal nulls are created. 

Phase D3 occurs due to the loss of one separator, which reduces the overlying 
field to multiplicity one, leading to four simply connected domains and only one 
separator. Thus, here T> = 4 and D — X = 3, in agreement with equation 1.1. In 
phase D4, the global double-separator bifurcation occurs twice. This creates four 
new separators bringing the total number of separators {X) to five. It also creates 
two source pairs of multiplicity three with two other source pairs of multiplicity 
one surrounding them. Here, therefore, V = J2n=i '"'^n = lx2-|-2x0-|-3x2 = 8 
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and once again equation 1.1 is satisfied. To create phase D5 the two source pairs 
of multipHcity three are reduced to multipKcity two each by the loss of a separator 
through the source plane. Thus, here V = X;^=i '^Dn = 1x2 + 2x2 = 6, X = 3 
and V — X = 3, so equation 1.1 is again true. 

6. Conclusions 

The magnetic skeleton is found for a magnetic configuration in which two oppo- 
site polarity photosphcric sources pass by each other in an overlying field. Such 
an elementary interaction represents a fundamental building block of the coronal 
heating process. Even though this magnetic field involves just four sources, two of 
which arc at infinity, the topology of the magnetic field is found to evolve through 
a series of relatively complex states with many multiply-connected source pairs and 
separators. The evolution found here in the dynamical MHD experiment is not 
dependent on the form of the resistivity used and has been found in experiments 
where 77 is uniform or of anomalous/hyper-resistive form and so is likely to be the 
natural evolution of the field. 

In general, it is expected that in magnetic fields with few sources finding the 
connectivity of the sources is the same as finding the fiux domains, i.e. it is expected 
that the volumes containing field lines of the same connectivity would be simple 
and have no holes. However, in our dynamic MHD experiment, we have found 
three phases (D2,D4 & D5) where one or more of the source pairs are not simply 
connected, but have more than one distinct route to link the same two sources. 
These multiply-connected source pairs arise naturally as a consequence of the global 
double-separator bifurcation. 

This new global double-separator bifurcation, which we have discovered here, 
is the way in which all the separators found in our dynamic MHD experiment are 
created. The bifurcation has two interesting consequences: 

• The new pair of separators engirdle a new coronal domain which is created 
when two separatrix surfaces touch and reconnection allows them to intersect 
above the source plane. This is likely to be one of the most natural ways in 
which separators are formed in the solar atmosphere. The new coronal domain 
pierces an existing domain making the source pair multiply-connected. Since 
multiply-connected source pairs are naturally created when a global double- 
separator bifurcation occurs they are likely to be common. 

• It appears as if the reconnection processes along each separator of a pair 
formed by the global double-separator bifurcation do not in general work to 
fill the same regions. Thus, for the newly formed domain to be created and 
to grow the reconnection along one separator must be dominant. 

In the potential evolution of the same magnetic setup, the magnetic field trans- 
forms through a non-generic transition state involving an isolated domain sur- 
rounded by an infinite number of separators. The dynamic MHD evolution has 
no such non-generic state, since the global double-separator bifurcation ensures 
that it is avoided. 

In the potential evolution there is generally only one separator. In the dynamical 
MHD experiment there are phases involving 1, 2, 3 or even 5 separators. The 
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reconnection is concentrated around these separators. Thus the separators are the 
focus for heating and so the detailed spatial distribution of the observed emission 
from this type of interaction could not be predicted by potential modelling nor, of 
course, could the temporal distribution of the thermal emission. 

The natural next questions arising from this work are, what is the nature of 
the reconnection and how rapid is the reconnection. The answer to these questions 
are not as straightforward as one might initially imagine and will be discussed in a 
follow-up paper (Parnell et al. 2007). 
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